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IC’  oa  <^oO  I ^The  results  of  ca)etilations  of  the  ionization  equilibrium  and  radiative  energy  ioaa  rates  Im 

C.  JI.  and  O iom  toufcw-denaity  plasmas  are  presented  for  electron  temperatures  In  the  range  <f0^ 
I ^ (gpT^  1^  eV).  The  ionization  structure  is  determined  using  the  steady-state  corona 

model,  in  which  electron  impact  ionization  ftom  the  ground  states  is  balanced  by  direct  radiative 
and  dielectronic  recombination.  Uaing  an  improved  theory,  detailed  calculations  are  carried  out  for 
the  dielectronic  recombination  rates  in  which  account  is  taken  of  all  radiative  and  autoionization 
prnrtiasfii  Involving  a single-electron  electric-dipole  transition  of  the  recombining  ion.  The  radiative 
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IONIZATION  EOUILIBRIUM  AND  RADIATIVE  ENERGY  LOSS  RATES  FOR 
C,  N,  AND  O IONS  IN  LOW-DENSITY  PLASMAS 


1.  INTRODUCTION 

An  important  method  of  determining  the  physical  properties  of  high-temperature  labora- 
tory and  astrophysical  plumas,  such  as  their  temperature  and  density,  is  the  analysis  of 
optically-thin  spectral  line  intensities  emitted  by  multiply-charged  atomic  ions.  In  addition 
the  presence  of  relatively  small  concentrations  of  multiply-charged  ions  in  a predominantly  hy- 
drogen plasma  can  have  an  important  effect  on  the  collisional  transport  processes,  such  as  parti- 
cle diffusion  and  thermal  conduction.  It  is  well-known  that  the  resonance  line  radiation  from 
incompletely  ionized  ions  can  influence  the  power  balance  in  laboratory  and  astrophysical  plas- 
mas as  well  as  provide  useful  diagnostic  information.  Accordingly,  spectroscopic 
investigations  have  been  conducted  to  determine  the  concentration  and  spatial  distribution 
of  multiply-charged  impurity  ions  in  such  plasmas. 

In  this  paper  we  present  the  results  of  calculations  of  the  distribution  of  ionization  stages 

and  radiative  energy  loss  rates  for  C,  N,  and  O ions  in  a low-density  plasma.  Only  a brief 

description  of  the  methods  will  be  given,  since  they  have  already  been  discussed  in  detail  in 

previously  reported  investigations  for  Fe  ions.  The  steady-state  optically-thin  corona 

model  is  applied  to  determine  the  distribution  of  ions  of  a given  element  among  the  various 

charge  states.  It  is  assumed  that  electron  impact  ionization  from  the  ground  states  is  balanced 

by  direct  radiative  and  dielectric  recombination.  The  ionization  equilibrium  distributions  thus 

obtained  are  independent  of  density  and  are  functions  of  only  the  local  electron  temperature. 

Since  the  low-lying  excited  states  which  are  produced  by  electron  impact  can  be  assumed  to 

undergo  radiative  decay  in  times  that  are  short  compared  with  the  time  between  collisions,  the 

emission  line  intensities  can  be  computed  from  the  corona  ionization  equilibrium  distributions 

of  (he  ground  sutes  and  the  ground  state  electron  impact  excitation  rates. 
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The  essential  innovation  in  the  present  investigation  is  the  use  of  an  improved  theory 
in  the  calculation  of  the  rates  for  dielectronic  recombination,  which  has  been  shown  to  be 
the  dominant  recombination  process  for  non-hydrogenic  ions  in  the  temperature  region  of 
maximum  equilibrium  abundance  and  maximum  resonance  line  emission.  The  simple,  and 
widely-used,  formula  of  Burgess  is  found  to  overestimate  the  dielectronic  recombination 
rates  for  certain  ionization  stages  because  of  the  neglect  of  the  effects  of  autoionization  into  an 
excited-state  of  the  recombining  ion. 

The  calculation  of  the  dielectronic  recombination  rates  is  discussed  in  Section  2.  In  Sec- 
tion 3 the  results  of  the  corona  ionization  equilibrium  calculations  are  presented.  The  theoreti- 
cal predictions  for  the  radiative  energy  loss  rates  are  given  in  Section  4. 

1.  DIELECTRONIC  RECOMBINATION 


Dielectronic  recombination  of  an  ion  A"*’*  with  residual  charge  z is  the  result  of  two  tran- 
sitions. The  first  is  a radiationless  capture  of  an  electron  accompanied  by  excitation  of  one  of 
the  bound  electrons  in  the  recombining  ion  to  form  a doubly-excited  state  / nl. 


(1) 

The  second  is  a stabilizing  radiative  transition  to  a single-excited  state  i,  nl  which  lies  below  the 
ionization  threshold  of  the  recombined  ion. 


JT  0.  «/)  - 

(2) 

The  dielectronic  recombination  rate  is  equal  to  the  product  of  the  capture  rate  and  the  branch- 
ing ratio  representing  the  probability  that  the  doubly-excited  state  decays  through  the  stabiliz- 
ing radiative  transition.  After  summation  over  all  important  excited  states  y of  the  recombining 
ion  and  quantum  numbers  nl  of  the  recombining  electron,  the  total  dielectronic  recombination 
rate  is  obtained  in  the  form 


2 


Rjiz.  l)  - N,  Nix.  i)  a^iz,  /). 


Jr 


(3) 

where  and  N(z,  i)  are  the  electron  and  initial  ((round)  state  ion  densities  and  a^iz.  i)  is 
the  dielectronic  recombination  rate  coeflTicient. 

The  branching  ratio  for  the  stabilizing  radiative  transition  is  obtained  by  dividing  the  ra- 
diative transition  rate  by  the  total  decay  rate  for  all  spontaneous  radiative  and  autoionization 
processes,  in  previous  calculations  for  the  dielectronic  recombination  rate  coefficients,  it  has 
been  assumed  that  autoionization  occurs  predominantly  into  the  ground  state  of  the  recombin- 
ing ion.  For  certain  excited  states  7 of  the  recombining  ion,  autoionization  into  a lower  excited 
state  can  occur  through  a S ft  •“  0 transition  with  a greater  probability  than  the  inverse  of  tran- 
sition (I).  The  inclusion  of  this  additional  autoionization  process  in  the  branching  ratio  for  the 
stabilizing  radiative  transition  produces  a substantial  reduction  of  the  dielectronic  recombina- 
tion rates  for  certain  ionization  stages. 

We  have  calculated  the  dielectronic  recombination  rate  coefficients  for  the  various  ions  of 
C,  M,  and  O,  taking  into  account  all  radiative  transitions  and  autoionization  processes  which  in- 
volve a single-electron  electric-dipole  transition  of  the  recombining  ion.  Because  of  the  fact 
that  large  values  of  the  outer-electron  principal  quantum  number  n play  the  most  imporunt 
role  in  the  dielectronic  recombination  of  low-z  ions,  the  autoionization  rates  can  be  obtained 
by  extrapolation  of  the  threshold  partial  wave  cross  sections  for  electron  impact  excitation  of 
the  recombining  ion.  With  increasing  density,  the  populations  of  the  most  important  nZ-levels 
will  be  significantly  influenced  by  electron  impacts’"’’  and  by  ion-produced  electric 
microfields.””  The  effective  recombination  rate  coefficient  defined  by  equation  (3)  will  be- 
come density-dependent.  The  investigation  of  this  density-dependence  will  be  the  subject  of  a 
forthcoming  paper. 
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3.  CORONA  IONIZATION  EQUILIBRIUM 


When  transport  processes  and  time- variations  can  be  neglected,  the  distribution  of  ions 
of  a given  element  among  the  various  charge  states  is  determined  by  the  steady-state  corona 
model  relationships 


Niz  - 1 )5(z  - 1 ) - N(z)a  (z),  1 < z < Z, 


(4) 


together  with  the  condition  that  the  sum  over  all  N{z)  correspond  to  the  total  abundance  N^. 
Here  S(z)  denotes  the  density  of  ions  in  the  ground  state.  5(z  — 1)  denotes  the  total 
rate  coefficient  for  electron  impact  ionization  from  the  ground  state  and  is  obtained  by  adding 
to  the  semi-empirical  direct  ionization  rate  coefficient  of  Lotz^'^^  the  contribution  describing 
autoionization  following  inner-shell  excitation.  Finally,  a (z)  is  the  sum  of  the  direct  radia- 
tive recombination  coefficient  and  the  dielectronic  recombination  coefficient  discussed  in 
the  preceeding  section.  Three-body  recombination  is  estimated  to  the  negligible  in  the  density 
region  of  interest. 


In  Table  1-3  the  relative  abundances  of  the  various  ions  of  C,  N,  and  Oare  presented  for 
electron  temperatures  in  the  range  4.0  < log|g  f,  CA)  < 7.0.  Similar  results  have  been  re- 
ported by  Jordan,  who  used  the  formula  of  Burgess  for  the  dielectronic  recombination 
rates.  The  inclusion  of  autoionization  into  an  excited  state  of  the  recombining  ion  reduces  the 
recombination  rates  for  the  lower  ionization  stages.  This  causes  a shift  of  the  maximum  abun- 
dance toward  lower  temperatures.  However,  this  shift  is  found  to  be  less  important  for  the 
low-z  ions  considered  in  the  present  work  than  for  the  high-z  Fe  ions  investigated 
previously.  This  is  a consequence  of  the  fact  that  the  2s  - 2p  transitions  of  the  recombin- 
ing ion,  which  are  not  affiscted  by  autoionization  into  an  excited  state,  are  found  to  be  most  im- 
portant in  the  low-z  ions.  The  present  results,  in  agreement  with  those  of  Jordan,  predict 
that  a given  ion  hu  its  maximum  equilibrium  abundance  at  a much  higher  temperature  than 
that  obtained  without  the  inclusion  of  dielectronic  recombination. 


The  contributions  to  the  collisional  transport  coefficient  from  multiply-charged  ions 
can  be  estimated  by  knowing  the  mean  charge  < z>  and  the  mean  square  charge  < > . 

The  mean  square  charge  also  enters  into  the  estimation  of  the  radiative  energy  loss  rate  due  to 
electron-ion  bremsstrahlung.  The  mean  and  mean  square  charges  obtained  for  C,  N,  and  O ions 
in  corona  ionization  equilibrium  are  shown  in  Figures  1-3  as  functions  of  temperature.  The 
< z>  and  < z^  > curves  show  initially  a rapid  increase  with  T,  which  is  followed  by  a flat 
region  where  the  //e-like  ion  is  predominant,  because  a substantial  increase  in  temperature  is 
required  for  ionization  of  a Is  electron  from  the  closed  shell.  The  effect  of  including  dielectron- 
ic  recombination  is  a substantial  reduction  in  the  values  of  < z > and  < z^  > at  a given  tem- 
perature. In  the  following  section,  it  will  become  apparent  that  this  shift  of  the  ionization- 
recombination  balance  in  favor  of  lower  ionization  stages  tends  to  enhance  the  line  radiation  in 
an  important  temperature  region. 

4.  RADIATIVE  ENERGY  LOSS 

The  power  radiated  per  unit  volume  as  a result  of  electron-ion  collisions  can  be  expressed 

as 


P «z. 


(5) 


where  the  efficiency  is  independent  of  density  and  is  a function  of  only  the  local  electron 
temperature.  In  corona  equilibrium  the  dominant  radiative  energy  loss  mechanism  is  expected 
to  be  resonance  line  radiation  excited  by  electron  impact.  The  efficiency  for  this  radiation  loss 
mechanism  is  given  by 


where  A£(z,  / 


‘z 

J)  are  the 


- Z L ^ 7). 

*-i  j (6) 

excitation  energies,  C(z,  i J)  are  the  collisional  excitation  rate 
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coefficients  obtained  from  a distorted  wave  calculation/'^^  and  S(z)IN2  are  the  ionization 
equilibrium  distributions  given  in  Tables  1-3.  We  have  taken  into  account  only  single-electron 
electric-dipole  transitions  from  the  valence  shell,  and  the  effects  of  self-absorption  have  been 

* 

neglected. 

During  the  process  of  dielectronic  recombination,  radiation  is  emitted  close  to  the  reso- 
nance line  of  the  recombining  ion.  The  efficiency  for  this  radiation  can  be  calculated  after  re- 
placing the  electron  impact  excitation  rate  coefficient  in  Equation  (6)  by  the  dielectronic 
recombination  coefficient  associated  with  the  i — j transition.  We  defer  to  a future  investiga- 
tion the  calculation  of  the  energy  loss  rate  for  the  radiation  which  is  emitted  when  the  n/-states 
of  the  outer-electron  cascade  to  the  ground  state  of  the  recombined  ion,  because  this  contribu- 
tion is  expected  to  have  a strong  density-dependence. 

The  efficiencies  for  electron  impact  excitation  of  resonance  line  radiation,  dielectronic 
recombination  radiation,  direct  recombination  radiation,  and  electron-ion  bremsstrahlung  are 
shown  in  Figures  4-6  for  C,  \,  and  O ions  Similar  calculations  have  been  carried  out  by  Cox  « 

and  Tucker,  but  the  individual  radiative  energy  loss  rates  were  not  presented. 

The  brood  maximum  near  10-  “K  (~  10  eV)  is  produced  mainly  by  2s  — 2p  transitions 
in  the  lower  ionization  stages,  while  the  smaller  peak  at  higher  temperatures  is  the  result  of 
Is  — 2p  transitions  of  the  //e-like  and  /Alike  ions. 

Although  the  radiation  emitted  through  the  stabilizing  transitions  in  the  dielectronic 
recombination  process  can  be  more  important  than  direct  recombination  radiation  and 
bremsstrahlung,  electron  impact  excitation  of  resonance  line  emission  is  clearly  the  dominant 
radiative  energy  loss  mechanism  for  incompletely  ionized  ions  On  the  other  hand,  the  shift  of 
the  ionization-recombination  balance  in  favor  of  lower  ionization  stages  which  results  from  the 
inclusion  of  the  dielectronic  recombination  rates  tends  to  enhance  the  loss  rate  due  to  line  ra- 
diation, which  falls  off  approximately  as  l/T,  at  the  higher  temperatures. 
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Table  1 — log|Q  N(z)IN^ioi  Cions 


log,oT,  CK) 

C I 

C II 

C III 

C IV 

c V 

C VI 

C VII 

4.0 

0.13 

0.56 

4.1 

0.86 

0.06 

4.2 

1.66 

0.01 

4.07 

4.3 

2.04 

0.00 

2.72 

4.4 

2.18 

0.01 

1.80 

4.5 

2.30 

0.03 

1.11 

4.6 

2.50 

0.12 

0.61 

4.54 

4.7 

2.82 

0.32 

0.29 

2.96 

4.8 

3.26 

0.61 

0.13 

1.77 

4.09 

4.9 

3.80 

0.99 

0.11 

0.90 

2.04 

5.0 

4.61 

1.63 

0.39 

0.46 

0.63 

5.1 

2.84 

1.27 

0.73 

0.11 

5.2 

4.22 

2.34 

1.29 

0.02 

5.3 

3.31 

1.81 

0.00 

5.4 

4.15 

2.26 

0.00 

5.5 

4.88 

2.64 

0.00 

4.00 

5.6 

2.97 

0.00 

2.61 

5.7 

3.27 

0.01 

1.50 

4.89 

5.8 

3.61 

O.IO 

0.68 

2.93 

5.9 

4.13 

0.40 

0.23 

1.56 

6.0 

4.86 

0.94 

0.15 

0.74 

6.1 

1.62 

0.33 

0.29 

62 

2.37 

0.66 

0.10 

6.3 

3.09 

1.03 

0.04 

6.4 

3.72 

1.38 

0.01 

6.5 

4.27 

1.68 

0.01 

6.6 

4.76 

1.95 

0.00 

6.7 

2.19 

0.00 

6.8 

2.40 

0.00 

6.9 

2.59 

0.00 

7.0 

2.76 

0.00 

I 


Table  2 - log,p  N(z)IN^[ox  Mens 


N I 

N 11 

4.0 

0.00 

2.42 

4.1 

0.05 

0.90 

4.2 

0.57 

0.13 

4.3 

1.43 

0.01 

4.4 

1.97 

O.OI 

4.5 

2.17 

0.02 

4.6 

2.32 

0.09 

4.7 

2.59 

0.26 

4.8 

2.99 

0.54 

4.9 

3.47 

0.89 

5.0 

4.04 

1.31 

5.1 

4.75 

1.84 

5.2 

2.59 

5.3 

3.85 

5.4 

5.5 

5.6 

5.7 

5.8 

5.9 
6.0 
6.1 
6.2 

6.3 

6.4 

6.5 

6.6 

6.7 

6.8 

6.9 
7,0 


N III  N IV  N V 


3.77 

2.32 

1.34 

0.71  4.46 

0.33  3.02 

0.15  1.97 

0.09  1.18  4.28 

0.14  0.62  2.63 

0.35  0.29  1.41 

0.80  0.28  0.64 

1.79  0.85  0.59 

3.14  1.84  1.04 

4.41  2.78  1.53 

3.61  1.96 

4.33  2.33 

4.96  2.65 

2.95 
3.30 
3.82 
4.52 


N VI  N VII  N VIII 


2.72 

1.03 

0.22 

0.04 

0.01 

0.00 

0.00 

3.57 

0.00 

2.32 

0.02 

1.33 

4.21 

0.12 

0.60 

2.50 

0.42 

0.23 

1.32 

0.94 

0.19 

0.61 

1.59 

0.38 

0.25 

2.28 

0.69 

0,10 

2.92 

1.02 

0,04 

3.50 

1.33 

0.02 

4.01 

1.60 

0.01 

4.47 

1.85 

0.00 

4.48 

2.06 

0.00 

2.26 

0,00 

9 


Ttble  3 - lof,g  N(i)IN ^ior  Olons 


CK) 

OI 

Oil 

OIII 

OIV' 

0 V 

0 VI 

0 VII 

0 VIII 

0 IX 

4.0 

0.02 

1.28 

4.1 

0.48 

0.17 

4.2 

1.56 

0.01 

7.50 

4.3 

2J9 

0.00 

5J0 

4.4 

2.80 

0.00 

3.45 

4.5 

2.99 

0.00 

2.20 

4.6 

3.14 

0.02 

1J2 

5.43 

4.7 

3J4 

0.09 

0.69 

3.64 

4.8 

3.66 

0.30 

OJO 

2J3 

4.9 

4.11 

0.63 

0.13 

1.41 

4.72 

5.0 

4.66 

1.06 

0.13 

0.76 

3.10 

5.1 

1.60 

0.27 

0J6 

1.89 

4.88 

52 

125 

0i8 

0.19 

1.03 

2.99 

5J 

3.03 

1.04 

0J4 

0.50 

1.60 

2.73 

5.4 

4.04 

1.76 

0.59 

0J4 

0.72 

1.02 

5.5 

3.02 

1.53 

0.84 

0.62 

0.23 

5 6 

4.61 

2.82 

1.75 

1.01 

0.05 

5.7 

4.04 

2.63 

1.45 

0.01 

5.8 

3.40 

1.85 

0.01 

4.22 

5.9 

4.08 

2.19 

0.00 

2.91 

6.0 

4.66 

2.48 

0.01 

1.88 

6.1 

2.73 

0.03 

1.10 

3.37 

6.2 

2.99 

0.13 

0.58 

2.02 

6.3 

3.31 

0.34 

0.33 

1.08 

6.4 

3.75 

0.70 

0.31 

0.49 

6.5 

4.32 

1.19 

0.49 

0.20 

6.6 

4.94 

1.75 

0.77 

0.08 

6.7 

2.28 

1.06 

0.04 

6J 

2.77 

1.34 

0.02 

6.9 

3.23 

1.58 

0.01 

7X) 

3.64 

1.80 

0.01 
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